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Abstract 

We show that the energy spectrum of the one-dimensional Dirac equation in the 
presence of a spatial confining point interaction exhibits a resonant behavior when 
one includes a weak electric field. After solving the Dirac equation in terms of 
parabolic cylinder functions and showing explicitly how the resonant behavior de¬ 
pends on the sign and strength of the electric field, we derive an approximate ex¬ 
pression for the value of the resonance energy in terms of the electric field and delta 
interaction strength. 
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1 Introduction 


Supercritical effects are perhaps one of the most interesting phenomena associated with 
the charged vacuum in the presence of strong electric helds [1,2] The study of super¬ 
critical effects induced by strong vector potentials goes back to the pioneering works of 
Pieper and Greiner [3], Zeldovich and Popov [4] among others. The idea behind super¬ 
criticality is to have positron emission induced by the presence of very strong attractive 
vector potentials. The phenomenon can be described as follows: the energy level of an 
unoccupied bound state dives into the negative energy continuum, i.e., an electron of 
the Dirac sea is trapped by the potential, leaving a positron that escapes to inhnity. 
The electric held responsible for supercriticality should be stronger than 2meC^, which 
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is the value of the gap between the negative and positive energy continua. Such strong 
electric helds could be produced in heavy-ion collisions [1,5]. A rigorous mathematical 
study of the behavior of the Dirac energy levels near the continuum spectrum and the 
problem of spontaneous pair creation has been carried out by Seba [6], Klaus [7] and 
Nenciu [8] among others. 

In order to get a deeper understanding of the mechanism responsible for supercriticality 
and for the resonant peaks appearing in the energy spectrum when supercritical helds are 
present, we proceed to work with a vector point interaction in the presence of a constant 
electric held. Point interactions potentials may be used to approximate, in a simple way, 
more complex short-ranged potentials. Among the advantages of working with conhning 
delta interactions we should mention that, they only possess a single bound state and 
the treatment of the interacting potential reduces to a boundary condition. The study 
of bound states of the relativistic wave equation in the presence of point interactions 
is a problem that has been carefully discussed in the literature [9,10,11,12,13]. The 
one-dimensional Dirac equation in the presence of a vector point delta interaction has 
also been a subject of study in the search of supercritical ehects induced by attractive 
potentials [1,14]. Soon after the publication of the paper by Loewe and Sanhueza [14], 
Nogami et al [15] pointed out that supercritical ehects are also absent in a class of 
non-local separable potentials in one dimension. 

Since we are interested in studying the mechanism of positron production by supercrit¬ 
ical helds, we proceed to analyze the resonant behavior of the energy when a bound 
state dives into the negative continuum. This resonant behavior is associated with the 
appearance of simple poles of the resolvent on the second sheet at a position very near 
the real axis [16]. 

The method of complex eigenvalues (Gamow vectors) was introduced in quantum me¬ 
chanics by Gamow [17] in connection with the theory of Alpha decay. Titchmarsch [18] 
and Barut [19] demonstrated an application, in the framework of non-relativistic quan¬ 
tum mechanics, for the Gamow vector method to the problem of an attractive delta 
interaction 6{r) of strength — cot a with a weak electric field term associated with the 
potential V (r) = —Ar. They found that the Schrodinger equation with a weak electric 
held exhibits a continuous spectrum from —cx) to -|-cx), and a resonance at E in the 
vicinity of Eq — lAtana, where Eq = — cot^ a is the energy of the unperturbed state. 
The eigenvalue E' is a complex number in the lower half plane. According to Barut, 
the Schrodinger equation with the potential —Ar describes a system that tries to form 
a bound state that “dissolves itself’ in the presence of the continuous spectrum. In this 
Letter, using the idea developed by Titchmarsch [18], we hnd the energy spectrum of 
the one-dimensional Dirac equation with boundary conditions associated with a vector 
Dirac delta interaction and a constant electric held whose strength is weak, and therefore 
it produces a perturbative ehect on the delta energy spectrum. We hnd that in this case 
the energy spectrum exhibits a resonance due to supercriticality. 

The Letter is structured as follows: In Section 2, we solve the one-dimensional Dirac 
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equation in the presence of an attractive S potential and a constant electric field. In Sec¬ 
tion 3, we compute the energy resonances and show how they depend on the electric field 
strength. We also derive an approximate analytic expression for the energy resonances. 
Finally, in Section 4 we summarize our conclusions. 


2 The one-dimensional Dirac equation 


In this section we will consider the 1-1-1 Dirac equation in the presence of the attractive 
vector point interaction potential represented by eV(x) = —g^S^x), and a constant 
electric held associated with the potential eV{x) = Xx. The Dirac equation, expressed 
in natural units {H = c = 1) can be written in the form [ 20 ] 

d' = 0, (1) 

where is the vector potential, e is the charge and m is the mass of the electron. 
The Dirac matrices 7 ^ satisfy the commutation relation { 7 ^, 7 ^} = with = 
diag(l,—1). Since we are working in 1 - 1-1 dimensions, we choose to work in a two- 
dimensional representation of the Dirac matrices 

7° = 0 - 3 , 7 ^ =(2) 

Substituting the representation matrix representation (2) into Eq. (1), and taking into 
account that the potential interaction does not depend on time, we obtain 


+ (Xx — E) + mas}X(x) 
ax 


= 0 , 


with T = a^X., and 

X(a;) = 

with the boundary conditions at x = 0 


^2, 


Xi(0+) = Wi(0 ) cos - 1X2(0 ) sinc/„, 
-^ 2 ( 0 ’^) = -iWi(0“) singf^ -h X2(0“) cosgy. 


( 3 ) 

( 4 ) 

( 5 ) 


The above conditions (5) describe a the point vector potential interaction of strength g^ 

[ 12 ]. 


Equation (3) is equivalent to the system of equations 

dX 

(m + Xx — E) Xi — i—r^ = 0, 

dx 


( 6 ) 
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dX 

+ (m — Xx + E) X 2 = 0, 
ax 

Introducing the new functions and ^2 

Xi = -|- iVl2i X 2 — III — iff 2 5 

we obtain that the system of equations (6) -(7) reduces to the form 

dVti 


dx 

dfl2 

dx 


+ i {Xx — E)Q,i — mD .2 = 0, 
— i {Xx — E)VL 2 — mVLi = 0, 


( 7 ) 


( 8 ) 


(9) 

( 10 ) 


which is more tractable in the search of exact solutions. Substituting (9) into (10) we 
obtain the second-order differential equation 


d^Ui 

dx"^ 


+ {iX + {Xx — eY — Ill = 0. 


( 11 ) 


Looking at the asymptotic behavior of the parabolic cylinder functions D^{z) [21] we 
obtain that the regular solutions, for A > 0, of Eq. (11) belonging to £^(—cx),0) and 
£^(0, cx)) {ImE > 0), respectively are 


A 


= {Xx-E) 

(x) 


( 12 ) 


BDp\\{-eX (A(r-E)l , 


where Dp y D_p_i are parabolic cylinder functions [21], p = and A and B are 

constants. 

Inserting (12) into (9) and using the recurrence relations for the parabolic cylinder 
functions [21], we obtain 


.V^ 


hlo (3^) = *- E'^AD_p \ —e *4 (Ax — E) 

m \ V A 




(13) 


From (12), (13) and (8), we have that the components of the spinor solution for x < 0 
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are 


Xf(x) = A 


X,-(x) = A 


D_p_i ( ^/-e *4 (Ax - E) 


'\/2A jZL p, I -ilL /< pN 

-e4D_p W-e 4 (Ax-E) 

m \ V A 


D_p_i I (Ax - i?) j + (Ax - i?) 


(14) 


analogously, we have that, for x > 0, the spinor X has the components: 


X+(x) = B 
Xt{x) = B 


Dp ( \l 4 (A X - E)j - 4 Dp_i 


'2 

-e*5 (Ax-£;) 
A 

-e*4 (Xx-E) 

A 


(15) 


Inserting (14) and (15) in (5), we obtain that the equation for the energy eigenvalues E 
has the form: 



The regular solutions of Eq. (11), for A < 0, can be obtained after interchanging the 
roles of r2b(^) r2^(x) and of (^) (^) Eq. (12) and Eq. (13) respectively; 

therefore Eq. (16) permits us to compute the energy resonances for positive as well as 
negative values of the electric held strength A 


3 Energy resonances 


It is not an easy task to solve Eq. (16), nevertheless it is possible to derive an analytic 
approximation for the energy eigenvalues. In order to do this approximation we consider 
the effect of the electric held in the vicinity of zero, i.e., where the delta interaction is 
located. We insert the spinor solutions obtained using this approach into (5), obtaining 
in this way an eigenvalue equation, whose roots are a good approximation to those of 
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(16). The approximate equation for fli has the form 


dx^ 


+ {i\ + E'^ — m^} fli = 0. 


(17) 


The solution of Eq. (17) that, for a; > 0, and A > 0 exhibit a regular asymptotic behavior 

is 


hi]*' = Aexp{—iViX + — m‘^x). (18) 

Substituting the solution (18) into Eq. (9) we obtain 

E exp{-iViX + -m^x), (19) 

m 

and, analogously, we have that for a; < 0, the corresponding solutions are 

rib = B exp{iViX + E‘^ — mXx), (20) 


n- = ^ m? + E ^ ^2 _ (21) 

m 

Substituting the right r2i(2 left r2j(2 spinor components into Eq. (5), we hnd that, 
for Qy = —TT, the energy state sinks into the negative energy continuum, exhibiting a 
resonant behavior that depends on the electric held strength as: 


where we have assumed that A is positive and small compared to m. It should be noticed 
that, when we turn off the electric perturbation, we recover the energy eigenvalue E(0) = 
—m for a point vector interaction with strength Qy = —tt [12]. It is worth mentioning the 
vanishing of the hrst order perturbative correction to the supercritical energy E = —m 
as a result of the presence of the electric held interaction A. This behavior can also 
be observed looking at the real part of the resonant energy given by Eq. (22). An 
approximate expression for the energy resonances for negative values of A can be obtained 
after substituting A by —A in Eq. (22). 

We can try to obtain an approximate solution to the full Eq. (11) after making a series 
expansion of the function iX + {Xx — E)^ — mX around zero. In this way, we have that 
the approximate function hli satishes the second order diherential equation 

u2q 

+ {iX +E"^-2XEx-m^}Vti = {). (23) 

The regular solutions fib for a; < 0, and for a; > 0 can be expressed, for A > 0, 
in terms of Hankel functions Edlj^^z) and [22] as follows 

nr = Az^'^Hh^iz), (24) 
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Fig. 1. Dependence of the real and imaginary parts of the energy on the electric field strength 
for 0 < A < 


Qt = (25) 

where A and B are constants and 2 : is given by the expression 


2; = 


(zA + E^-m^- 2E\xf^ 
3Xe 


(26) 


Substituting the solutions (24) and (25) for 12^ and fl]*' into Eq. (9) we obtain the 
corresponding expressions for and Taking into account the relation between 
Xi, X 2 and rii, (8), and inserting the lower and upper components of X into the 
boundary condition (5), we find that an approximate equation for the energy spectrum 
is: 


(a-2(0-)A-,(0+) - A2(0+)Ai( 0-)) cosg„ + i (A2(0+)A2(0+) - A,(0+)A,(0-)) sing„ = 0 

(27) 

It is worth noting that the energy spectrum obtained from Eq. (27) with and ni*" 
given by Eqs. (24) and (25) gives a result comparable to the one obtained using the 
relation Eq. (22) Fig. 1 shows how the real and imaginary part of the energy depend on 
the parameter A when a point vector interaction of strength = —tt is applied. The 
energy exhibits a resonant behavior produced by the electric held. 

The inverse of the imaginary part of Eq. (22) permits one to determine the time [23,24] 
that the state stays in negative continuum. The stronger the electric held perturbation, 
the shorter the mean life of the resonance [23]. An interesting result that can be observed 
from Eq. (22) and it is also depicted in Fig. 1 is that the imaginary part of the resonant 
energy behaves as a square root with respect to the real part of the energy. 
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4 Concluding remarks 


In this article we have shown that the presence of a weak electric held plays a crucial 
role in the appearance of a resonant energy state in the one-dimensional Dirac equation 
in the held of an attractive vector delta interaction. This phenomenon is analogous to 
the Stark ehect [25], where a perturbative constant electric held produces a shifting and 
broadening of the energy levels of the hydrogen atom. It also parallels the autoionization 
process, where the presence electron-electron interaction term induces a self-ionization 
process and levels decay like in the Auger ehect [16,25]. Such analogies seem to indicate 
that, for resonant ehects, the shape of the perturbative potential is sometimes more 
important than the strength of the potential held. 

We have also shown that, the one-dimensional Dirac equation in the held of an attractive 
vector delta interaction exhibits a resonant supercritical behavior when one introduces 
a potential of the form Xx. Since this ehect is not present when one consider the sole 
vector interaction —g^5{x) [12], we conclude that it is necessary the presence of the 
electric held in order to sink the delta bound state into the negative continuum. 
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